The effect of a radome 
on submillimeter site testing measurements 

Paolo G. Calisse 

School of Physics, University of New South Wales, Sydney NSW 2052, Australia 
e-mail: pcalisse@phys.unsw.edu.au 

KEYWORDS 

atmospheric effects, site testing, submillimeter 

ABSTRACT 

We evaluate the effect that radome transparency has on atmospheric 
opacity measurements performed by the skydip technique. We show that, 
except at rather high opacities, it is not sufficient to ignore losses in the 
radome (or "window") during the data analysis and then subtract them off 
from the derived atmospheric opacity. Perhaps surprisingly, unless radome 
transparency is correctly modelled, the atmosphere will appear to have a 
minimum opacity that is many times greater that the radome losses. Our 
conclusion is that some previous site studies may have significantly 
underestimated the quality of the best submillimeter sites, and that the 
difference between these sites and poorer sites may be much greater that 
currently believed. We also show that most of the residual 857 GHz 
opacity at the best sites, currently ascribed to "dry-air opacity", can in fact 
be just an artefact caused by not properly modelling the radome during the 
data analysis. 

INTRODUCTION 

The skydip technique (Dicke et al. 1946) is used to measure the atmospheric opacity from 
radio to infrared wavelengths, both for site testing purposes and to calibrate observations. 
The technique requires a measurement of the sky flux at different elevations, and derives 
the zenith sky opacity on the assumption that the atmosphere can be modelled as a single 
slab at uniform temperature. 

Skydips have been widely used in the sub -millimeter and millimeter wavelength region 
(Chamberlin & Bally 1994; Chamberlin, Lane, & Stark 1997; Dragovan et al. 1990; Lane 
1998; Matsuo, Sakamoto, & Matsushita 1998; Pardo, Serabyn, & Cernicharo 2001; 
Radford & Holdaway 1998), because, when observing astronomical sources at these 
wavelengths, the attenuation introduced by the atmosphere, and its fluctuations, are 
amongst the main sources of noise, calibration error, and systematic errors. 
We have reviewed the different algorithms that have been used for data reduction, which 
allow for different contributions. In the general case one should include terms related to 
the sky, such as the cosmic background radiation, or to the telescope and receiver, like 
spillover, blockage and ohmic losses (Ulich & Haas 1976). With a finite beam size the 
non-linearity of the relation between the elevation and the airmass number when 
approaching the horizon should also be considered. This limits the minimum elevation to 
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use for the skydip depending on the instrument beam size (Rohlfs & Wilson 1996) and 
the geographical characteristics at a particular site. 

However, the algorithms traditionally used do not always take into account the effect of 
absorbing material used for the radome or window 1 . We will show that this effect cannot, 
in many conditions, be treated as a simple offset to be subtracted later from the sky 
opacity. Moreover, in particularly good atmospheric conditions the instrument will be 
unable to obtain a valid and unambiguous value for the sky opacity. 
This ambiguity can, however, be fully removed using an appropriate data reduction 
algorithm, as shown below. 

THE ATMOSPHERIC SINGLE SLAB MODEL 

Most authors make use of the standard atmospheric model (in which the atmosphere is 
considered to be a plane -parallel slab at uniform temperature) to describe the relation 
between the atmospheric emission and the airmass number X = l/sin(z) at elevation z. 
Assuming that one is using a chopping technique, the signal at the detector (expressed in 
temperature units), r'(^0,is: 

T'(X) = T' alm (l-e- xr »)-T ref (1.1) 

where t' is the estimated zenith opacity, T' alm the physical temperature of the 
atmosphere, and T ref the reference temperature against which the chopping is performed. 

All the other corrections (cosmic background radiation, spillover, ohmic losses, etc.) have 
been neglected. A prime will be used to distinguish values obtained neglecting the 
window transparency from those that properly include this correction. 

OPACITY OBTAINED FROM MEASUREMENTS 
AT 1 AND 2 AIRMASSES 

From eq. (1.1), we obtain at 1 and 2 airmasses: 

T'(l) = T' atm (l-e-*°)-T ref (2.1) 

T'(2) = T' al Jl-e- 2t '°)-T ref (2.2) 
We can analytically solve this system of two equations for % 1 and T \ tm , as follows: 
, r(2)-r(l) 



T'(V + T ref 
T'(D + T ref 



(2.3) 



T' alm =— ^ (2.4) 

1-y 

where y ' = e~ x " is the atmospheric transmission. 



1 In some cases the effect of the radome can be correctly modeled or measured as a component of the 
antenna beam coupling efficiency. 
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However, if a radome is used, the detected signal T(X) at a given airmass number X 
should properly be written as: 

T{X) = » win T am (1 - ) + T win (1 - » win ) - T ref (2.5) 

where & win is the transparency of the window, and T win its temperature (De Zafra 1995). 

We have neglected scattering and reflection by the window for the following reasons. 
Scattering would essentially introduce a spillover term affecting the relation between the 
model and experimental data, especially at low elevation, where the airmass number 
changes dramatically for little change of elevation. This deviation is not seen in real data 
(Calisse et al. 2002). Reflection would introduce a thermal signal in the same way that 
absorption does, and can be treated simply as another form of absorption as long as all 
temperatures in the vicinity of the window are comparable. We show below that the 
window temperature is in fact not critical, and hence this is a good approximation. 

For the case of Gore-tex, a fiber widely used as a radome material, these assumptions 
are also supported by direct measurements. 
Within this approximation, we have at 1 and 2 airmasses: 

m = ^J atm (1 - e- z " ) + T win (1 - V win ) - T ref (2.6) 
7X2) = V w J am (1 - e 1 ^ ) + T mn (1 - & win ) - T ref (2.7) 

The correct value of the zenith opacity to can be found substituting the values of T at m 
obtained from equation (2.6) into eq. (2.7): 

y = r(2) - r(1) (2.8) 

T(l)-T win -(l-# win ) + T ref 

T _ T(D-T wm (l-i} wm ) + T ref 

tf (l-y) K } 

where y = e~ T ° is now the actual atmospheric transmission. 

We will now evaluate the error that is introduced by ignoring radome losses, as in eq. 
(1.1), substituting the values for T(l)and T{2) obtained from eqs. (2.6) and (2.7) into the 
approximate solutions (2.3) and (2.4). 
We obtain: 

y' = ^y (2.10) 

k-y 

(^ )^ 

T atm ~ ^nirr^atm ~, Z ' 2 (2 - 1 1 ) 

k-2y + y 
T 1-1} 

where k = 1 + — . The apparent zenith opacity T ' is plotted against the actual 

opacity T , for typical conditions, in Fig. 1. 
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In high opacity conditions the difference between the actual and the apparent opacity is 
independent of the opacity value, as expected. In fact, for T » 1 , we can neglect y with 

respect to 1 in eq. (2.10) and obtain an asymptotic correction for T and T am : 

T =T' -ln[£] (2.12) 

^i»ji = TT T atm (2.13) 

Note that k > 1 ; that is, neglecting the window opacity always leads to an over-estimate 
of the atmospheric opacity. 

For & win = 0.81 , a typical Gore-Tex transmission at 857 GHz, (Radford 2002), 
T win = 210K and T atm = 22QK (typical winter South Pole conditions) we 
haveT =t' -0.25 and T atm =0.96-T' alm for high opacity values. 

However, in the more important case of medium and low opacities, the approximate 
model fails altogether. In fact, there is a minimum value to the apparent opacity that the 
skydip technique will measure when using a window of transparency & win . This value can 
be obtained by finding the zeroes of the derivative of y' with respect to y: 

y mm =k-Jk-(k-l) (2.14) 

For a given radome transparency, as the sky opacity progressively decreases, the 
approximate model creates first an increased offset, then a reduced sensitivity to opacity 
changes, and, finally, the apparent opacity actually increases, and the model described by 
eq. (1.1) generates an ambiguity in the opacity data. 

If the correct model is used to solve the system of equations or to fit the data, this 
problem disappears. 

To understand when, in general, the ambiguity can occur, we plotted the relation between 
the radome transparency and the apparent minimum opacity (see Fig. 2), for two 
different ratios between the radome temperature T w t n and the actual atmospheric 
temperature T at m. We find that the results are not very sensitive to the ratio between these 
two temperatures. This can be explained by the fact that the emission term in eq. (2.5) is 
usually negligible with respect to the absorption term. So, even in the presence of 
relatively high (30-50 K) uncertainties in the window temperature we will still obtain 
valid results. 



OPACITY OBTAINED BY NON-LINEAR FITTING 

In the previous section we discuss measurements taken only at 1 and 2 airmasses. 
However, opacity measurements are often made by observing at many different airmasses 
and using a non-linear fit to the data. As demonstrated in another paper (Calisse, et al. 
2002), we observed an excellent agreement between the parameters obtained with both 
techniques for large sets of opacity data obtained at Dome C and South Pole. 
For data already published in the literature that did not properly account for the radome 
opacity, we can nevertheless calculate a correction. The correction to be applied will be 
now slightly sensitive to the details of the sampling (number and values of airmass steps) 
and so obtaining an analytical expression for the correction will not be possible. 
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For this reason, to correct published data we simulated the apparent atmospheric 
brightness temperature at different airmasses. Then, we applied a non-linear fitting 
obtaining the apparent zenith opacity and atmospheric temperature. 
The results are not significantly different from the 1-2 airmass analytical solution: a 
minimum apparent opacity has still been found. See Fig. 3, where we simulated an actual 
observation where skydips were made with measurements at seven different airmasses 
from 1 to 4 and different opacity values ranging from 0.1 to 4. 

We obtained a minimum apparent opacity of 0.6, corresponding to an actual opacity of 
about 0.2, in good agreement with the minimum value reported at 857 GHz in sites like 
the South Pole (Radford 2002). The correction to be introduced at high opacity values has 
also been determined to be about 0.3. 

Note that as the response of the instrument is not very steep around y min , a small 

variation in the apparent opacity can correspond to a large variation in the actual value. 
For this reason, statistics on the atmospheric noise obtained in these conditions should be 
treated with caution. 

Moreover, the value of the reduced^ 2 value rises to very high values at low opacity (for 
example, 24.5 for T = 0.1 in the simulation described above), indicating the poor 
capability of the model (1.1) to explain the observed data and noise distribution. 

CONCLUSIONS 

A non-linear error, which becomes very serious at low opacity, is introduced by data 
reduction algorithms that do not properly take into account the radome transparency. 
The consequences are as follow: 

• At high opacity values, even if the original data are no longer available and the 
approximate model (1.1) has been used in data analysis, it is still possible to 
correct the data by making use of the correcting algorithm (2.10) and (2.1 1). 

• The exact model should be used whenever the apparent atmospheric opacity is 
equal to or lower than the one given by Fig. 2 at the particular window 
transparency, in order to remove a possible ambiguity. For example, at 857 GHz 
(T at m=220 K, Twi„= 270 K, transmission=0.81) the exact model (2.5) should be 
used whenever the apparent opacity at a site falls below ~1. In fact, in this case 
the approximate model (1.1) fails to produce correct data and statistics. Corrected 
results at three relevant sites (Radford 2002) are shown in Table 1. Similar 
corrections apply to other submillimeter windows. 

• While the extent of the radome effect depends strictly on the transmission 
properties of a specific radome material, it is not very sensitive to the radome 
temperature, as the emitting term is negligible with respect to the absorbing one in 
eq.(2.5). 

• This "radome effect" can explain the "hard lower limit" to the apparent sky 
opacity observed on some skydip data series obtained at South Pole and at other 
sites with particularly good atmospheric conditions (Radford 2002; Radford & 
Holdaway 1998). Significantly, the minimum detectable opacity of 0.6 derived 
from Fig. 2 corresponds to the apparent minimum detected opacity for three sites 
(Chajnantor, South Pole, Mauna Kea). 



5 



• The statistics at sites where the opacity is very low (i.e. South Pole in winter) are 
affected more than those at other sites. For this reason, South Pole could be 
relatively better than previously reported in terms of average opacity at 
submillimeter wavelengths. 

• Noise statistics would be affected by the window transmission, but in the opposite 
way: a change in the real opacity does not induce significant variations in the 
apparent one, so that atmospheric opacity variations will be underestimated. 
However, this effect is partially mitigated, as seen in statistics on actual data 
obtained at Dome C (Calisse, et al. 2002), by the fact that the correct model (2.5) 
fits the real data better, reducing the uncertainty and, in turn, data dispersion. 

• Several relations between opacity at submillimeter wavelengths and PWV, 
obtained by correlating the values of PWV obtained by onsite radiosonde 
measurements and opacity apparent by radiometers (Chamberlin 2001) have been 
found. In general, a consistent residual opacity at very low PWV values has been 
observed and ascribed to a dry-air component. However atmospheric models 
(Pardo, et al. 2001) predict a much lower level for the dry-air opacity than that 
measured, particularly at 857 GHz. We believe that this "dry" component has 
been greatly overestimated in some cases by not correctly taking into account the 
window properties. 
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Fig. 1 - Apparent atmospheric opacity as a function of the actual opacity for a range of 
different window transparencies, when using the analytical solution available for 1 and 2 
airmass measurements. In this simulation T. = 220^ and T. = 270^ . 
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Fig. 2 - Minimum detectable opacity, obtained by the analytical solution of measurements 
at 1 and 2 airmasses, as a function of the radome transparency. The continuous line 
represents the case of T win lT atm = 1 , the dashed line represents the case of T win lT alm = 1.5 . 
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Fig. 3 - Apparent opacity as a function of actual opacity for the case where a non-linear 
fitting of data to the uniform slab atmospheric model has been used, but the window 
transparency ignored. In this simulation is T atm = 220K and T win = 210K , as in the case of 
1-2 airmass (see Fig. 1). 
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Broadband Zenith Opacity at 350 um 
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Chajnantor 


1.21 


0.95 


1.68 


1.4 


2.52 


2.3 


Mauna Kea 


1.46 


1.2 


2.15 


1.9 


3.14 


2.9 


South Pole 


1.23 


1.0 


1.52 


1.3 


1.94 


1.7 



Table 1 - Recalibration of the cumulative percentile distribution of broadband opacity at 
three relevant sites for submillimeter astronomy, by applying the correction algorithm to 
results found in the literature (Radford 2002). x \ is the published zenith opacity and T the 
actual zenith opacity. The parameters used calculating the correction are shown in Fig. 3. 
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